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Theoretical attenuation and phase characteristics at VLF are presented for a number of 
idealized models of the lower ionosphere. The results indicate the limitations of the sharply 
bounded model. 



1. Introduction 

Despite the great number of papers dealing with 
VLF radio propagation, there appears to be a need 
to consolidate some of the recent theoretical data. 
It is the purpose of this paper to fulfill this need. 
The intent is to present results, in convenient 
graphical form, which illustrate the expected value of 
the attenuation and phase velocity of the dominant 
modes in the frequency range from 8 to 30 kc/s. 
Some of the effects considered are due to the earth 
curvature, ground conductivity, and finally the 
gradient of the effective conductivity of the lower 
ionosphere. 

No attempt will be made here to describe the 
mathematical aspects of the subject in any detail. 
However, certain basic formulas will be introduced 
to permit comparisons with previous work [Wait, 
1962]. 

The basic theoretical model is a smooth spherical 
earth of homogeneous conductivity surrounded by a 
concentric (spherically) stratified ionosphere. Ini- 
tially the ionosphere is regarded as a sharply bounded 
isotropic medium. For this discussion the influence 
of the terrestrial magnetic field is not considered. 
The non-sharpness or diffuseness of the lower edge 
of the ionosphere is treated by using an exponential 
isotropic model which has the virtue of simplicity. 

2. Basic Modal Equation 

To permit a straightforward treatment, it is as- 
sumed that the level h above the ground may be 
characterized by an ionospheric reflection coefficient 
Rj. Furthermore, it is convenient to make the 
substitution 

R t =-ex V (olC), (1) 

j where 0" is the cosine of the angle of incidence at the 
reference level in the ionosphere. While the a is 
generally a (complex) function of C", it is found that 
under most practical situations a is approximately 
a constant for C ranging from up to about 0.3. 



This, of course, encompasses the important angles in 
long distance propagation of VLF radio waves. 

The relevant modal equation for the low order (im- 
portant) TM modes may now be written [Wait, 1962] 
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and 



&=27r/wavelength, 
a=radius of the earth, 
n=mode number, 

t=-(ka/2) 2/ *C\ 



This form of the mode equation is valid for a ground 
surface which is effectively a perfect conductor. The 
logarithmic term involves the derivatives of the 
Airy functions of argument t. 

Under the somewhat unrealistic assumptions that 
the ionosphere is a sharply bounded ionized medium 
and the terrestrial magnetic field may be neglected, 
it is known that [Wait, 1962] 
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where oi r =^ 2 jv. Here 

co =the (angular) plasma frequency, 
y=the effective collision frequency. 



(3) 



by 



The attenuation A n , of the nth. mode, is then given 



A n =— Im kS n nepers/meter, 



where #»=(1 — C2)*^l — CJ/2. In more convenient 
units 
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Figuee la. Attenuation rate of the first mode for a perferctly 
conducting earth and a sharply bounded ionosphere. 
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Figure lb. Attenuation rate of the second mode for a perfectly 
conducting earth and a sharply bounded ionosphere. 



A n ^lm kC 2 J2XS.mxiO z in db/1000 km, (4) 



when 



ft=2^ kc/B /300, 



where / kc/s is the frequency in kilocycles per second. 
The phase velocity v n of the nth mode, referred 
to the velocity of light in vacuo, is then obtained 
from 



^-l^Ke(C 2 J2). 



(5) 



The complex values of C n are the roots of the modal 
equation quoted above. It is considered to be con- 
venient to number the modes in the VLF range so 
that mode number 1 is the mode of least attenuation. 
Higher order modes correspond to modes of progres- 
sively greater attenuation. 

3. Attentuation for the Sharply Bounded 
Model 

Using (2) and (4), the attenuation rates A x and 
A 2 for the first two modes of lowest attenuation are 
shown in figures la and lb, respectively. The height 
of the reflecting layer varies from 60 to 100 km in 
5 or 10 km intervals. The ground conductivity a g is 
here taken to be oo while co r has the value 2 XIO. 5 
The latter is typical of daytime conditions. It is in- 
teresting to note that the attenuation A\ for the first 
mode has a broad minimum when plotted as a func- 
tion of frequency. It is also noted that at the 
lower frequencies the attenuation A x increases rather 
significantly for the lower reflecting heights. 

In general it is indicated that the attenuation A 2 
of the second mode is somewhat greater than A 1 . 
It is important to observe that A 2 is very dependent 
on the height parameter h. 

The influence of finite ground conductivity can be 
determined by solving a slightly more complicated 
form of the modal equation. To within a good 
approximation it has the same form as (2) if the log 
term is replaced by 
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where 



q=-i(ka/2)lZ g /120T, 



and where Z g is the surface impedance of the ground. 
Choosing the ground conductivity a g to have a 
number of typical values, A x was calculated from 
the modified form of the mode equation. It appears 
that for propagation over sea water (i.e., a-^4000 
millimhos/m) the value of <r g is effectively infinite. 
However, over moderately or poorly conducting 
land, the effect of finite ground conductivity cannot 
be ignored as indicated clearly in figure lc. 
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Figure lc. Attenuation rate of the first mode for an imperfectly 
conducting earth and a sharply bounded ionosphere. 
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Figure 2a. Phase velocity of the first mode for a perfectly 
conducting earth and a sharply bounded ionosphere. 

4. Phase Velocity for the Sharply Bounded 
Model 

The values of the phase velocity of the modes 1 
and 2 for the same conditions as in section 3 are 
shown in figures 2a, 2b, and 2c. As expected, the 
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Figure 2b. Phase velocity of the second mode for a perfectly 
conducting earth and a sharply bounded ionosphere. 
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Figure 2c. Phase velocity of the first mode for an imperfectly 
conducting earth and a sharply bounded ionosphere. 



phase velocities differ only slightly from c, the veloc- 
ity of light in vacuo. Nevertheless, this difference 
is very significant. It is interesting to note that the 
phase velocity Vi of the first mode may be actually 
equal to c for frequencies around 15 kc/s when the 
reflecting heights are around 65 km. For frequen- 
cies above this, the phase velocity is less than c 
and the mode becomes a "slow wave." As indicated 
by figure 2c, the ground conductivity appears to 
have only a very small influence on the phase 
velocity of VLF radio waves. 
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5. A Simplified Form of the Modal Equation 

For certain extensions and other applications, it 
is important to ascertain if any simplifications can 
be made to (2). The key step is to recognize that 
when t is near zero, the factor i log [iv^t) /w[(f)] 
may be replaced by 7r/3. This is valid when \{ka/ 
2)* C\ is small compared with unity. At the same 
time, terms such as (C 2J r2h/a)* are replaced by 



(2A/a)*(l- 



C 2 a\ 



which is justified when C 2 is somewhat less than 
h/a. Employing these approximations, the modal 
equation (2) reduces to a simple algebraic equation 
in C 2 . Explicitly, the solution may be written 



/12 

^ n 



1271-5] (7r/6)-(2&a/3) (2h/a) s/2 -ia(2h/a^ 



ka(2h/a)> + (ia/2) (2h/a)~ 



(6) 



To give some idea of the goodness of this simple 
formula, the phase velocity i\ for the first mode is 
shown plotted in figure 3 for a=0. It is quite 
evident, that for most purposes, the approximate 
method of solution is adequate. 



>-*■ 

o 



< 
x 
a. 



0,006 



0.004 



0.002 



0.000 



-0.002 



-0.004 



-0.006 



-0.008 



-0.010 



COMPARISON OF EXACT (— 

AND APPROXIMATE ( ) 

METHODS OF SOLUTION 




16 



20 



28 



FREQUENCY, kc/s 



Figuke 3. Comparison between accurate or "exact" modal and 
an approximate form used for simplified calculations. 



6. Exponential Layer 

To illustrate the influence of a non-sharp bound- 
ary at the lower edge of the ionosphere, it is assumed 
that the parameter a> r varies exponentially with 
height z. Thus 



oo r =0Q r exp 



((8(*-A)) f 



(7) 



where co r is the value of oo r at the reference height 
z=h. For a t}^pical D-layer profile, the exponen- 
tial form is a surprising^ good fit if /3 is about 0.5 
km -1 . For daytime, h is approximately 70 km while, 
for nighttime conditions, h is about 90 km. For 
these calculations u r is assigned the value 2.5 X10 5 . 

The reflection coefficient for such an exponential 
layer may also be expressed in the form — exp 
(aC) when referred to the reference height z=h. 
Furthermore, a is again nearly a constant for a wide 
range of angles of incidence. A numerical treat- 
ment of the exponential layer for plane wave in- 
cidence has already been given in some detail 
[Wait and Walters, 1963]. For present purposes, 
the complex values of a previously obtained are used 
in the simplified modal equation (6) given above. 

For presentation of results, the attenuation and. 
phase velocity are presented as differences to the 
corresponding results for the sharpfy bounded model. 
This technique has two advantages: (1) it makes 
use of the more accurate attenuation and phase 
data using the sharply bounded model, and (2) it 
illustrates the influence of a non-sharp boundary in 
a direct manner. 

The attenuation increment A A n , which is plotted 
in figures 4a, b, c, and d, is defined by 



AA n =A n W-A n , 



(8) 



where A n { P is the resultant attenuation rate for the 
n th mode when the ionosphere has an exponential 
profile, while A n is the attenuation rate of the sharply 
bounded model with co r = 2X10 5 . 

The attenuation increment A A x for the first mode 
is shown in figures 4a and b for h values of 70 and 
90 km. The corresponding attenuation increments 
A A 2 for mode 2 are shown in figures 4c and d. It 
is interesting to observe that, for the larger values 
of 13 corresponding to a relatively abrupt boundary, 
AA n is relatively small and does not depend signif- 
icantly on frequency. For smaller values of ft, 
corresponding to a diffuse boundary, the attenuation 
rate may be significantly increased. This is particu- 
larly apparent at the higher frequencies. 

The phase velocity increment N is defined in a 
similar fashion. Explicitly, 

N n =\~% (9), 

where v^ ] is the phase velocity of the nth. mode 
for an exponential ionosphere while v n is the cor- 
responding phase velocity for the sharply bounded 
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Figures 4a, b, c, and d. The increase or increment of the attenuation for a diffuse ionosphere over that of a sharply bounded 

ionosphere. 

(The effective conductivity of the diffuse ionosphere varies according to exp [fi(z-h)].) 
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Figures 5a, b, c, and d. The increase or increment of the phase velocity for a diffuse ionosphere over that of a sharply bounded 

ionosphere. 
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model as indicated in figure 5. It is apparent that 
N n is a small quantity for the range of frequencies 
indicated. In general, it may be observed that a 
diffuse ionosphere (i.e., small /3 values) is associated 
with higher values of phase velocity. Physically, 
this is a consequence of the lowering of the effective 
height of reflection. It may be noted that the effect 
is particularly noticeable for the higher range of 
(VLF) frequencies. 

The actual ionosphere is expected to encompass 
the whole range of fi values indicated in figures 4 
and 5. Thus, the phase perturbations are a sensitive 
indicator of the relative diffuseness of the lower 
ionosphere. 

The reader should be cautioned that the results 
in the present paper make no allowance for the 
terrestrial magnetic field. This factor will change 
the attenuation rates to some extent but the con- 
clusions regarding the gradualness of the lower 
ionosphere are still valid. 
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